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SUMMARY 

A procedure for the local stiffness modifications of large structures is described. It enables 
structural modifications without an a priori definition of the changes in the original structure and 
without loss of efficiency due to multiple loading conditions. The solution procedure, implemented 
in NASTRAN, involves the decomposed stiffness matrix and the displacement vectors of the original 
structure. It solves the modified structure exactly, irrespective of the magnitude of the stiffness 
changes. In order to investigate the efficiency of the present procedure and to test its applicability 
within a design environment, several real and large structures were solved. The results of the effi- 
ciency studies indicate that the break-even point of the procedure varies between 8% and 60% 
stiffness modifications, depending upon the structure’s characteristics and the options employed. 

INTRODUCTION 

For the efficient design of large and complex structures it is necessary to investigate several 
modified structural configurations. Frequently, structural modifications are made at a few discrete 
locations in the structure, affecting only a relatively small part of it. Consequently, it is desirable 
to apply a modification procedure which enables rapid re-analysis of the structure in a time step 
which is considerably less than that of the analysis of the regular structure. Furthermore, the pro- 
cedure should be automated with minimum requirements of input data. 

The various modifications procedures which have been developed so far can be subdivided 
into iterative and direct procedures. The iterative procedures are based either on a simple technique 
of successive approximations (reference 1 ), or on a version derived from the Gauss - Seidel technique 
(reference 2), or on a first order Taylor series expansion of design variables (references 3 and 4). 

The procedures become less efficient as more loading conditions are applied to the structure and 
more inaccurate as the magnitude of the stiffness changes increases. The advantage of using these 
procedures is that the efficiency of the procedure is not affected by the relative number of modi- 
fied degrees of freedom. The direct procedures (references 1 and 5), are exact and their efficiency 
is mainly governed by the relative number of the modified degrees of freedom. Their efficiency is 
not affected by the number of loading conditions, and the accuracy of the analysis results is inde- 
pendent of the magnitude of the stiffness changes. Argyris and Roy (reference 6) have developed 
a general direct modification procedure which handles coupled combinations of three types of 
modifications, namely, changing stiffnesses of elements, adding degrees of freedoms and removing 
degrees of freedom. 
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To our knowledge, the direct procedures have not yet been investigated for large and com- 
plex structures and their applicability within a design environment has not been tested. In addition, 
operation counts alone, as it was done in references 1 , 5 and 6, are not sufficient to predict the 
overall efficiency of the procedures. This is because, for large and real structures, the speed of in- 
core operations relative to the speed of out-of-core data transfer, the type of finite elements and 

the time required for the assembly of their stiffness matrices may have a considerable effect on the 
efficiency. 

The present paper describes a direct procedure for the solution of structures with local 
stiffness modifications (reference 7). The basic concepts are similar to those of Argyris and Roy 
(reference 6). However, their procedure was modified and implemented in the static solution of 
NASTRAN (reference 8) by employing NASTRAN'S DMAP language (references 9 and 10). In 
order to reduce the burden of input data preparation, a NASTRAN preprocessor program, REAN, 
has been developed. This preprocessor, in addition to its other features, takes care of the charac-’ 
teristics of the modified and unmodified structure constructing the Boolean transformation matrices, 
which relate their corresponding degrees of freedom. 

In order to investigate the efficiency of the modification procedure and to test its applica- 
bility within a design environment several real and large structures were solved. For each of the 
examples a number of modifications were made, including the addition or removal of elements 
and the investigation of regions with reduced stiffnesses in order to study the post-buckling behavior 
of these regions. The size of the different mathematical models varied between 60 and 1 680 un- 
constrained degrees of freedom. 
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SYMBOLS 

Boolean transformation matrix 

inverse of [K] 

modulus of rigidity 

unit matrix 

stiffness matrix 

modified stiffness matrix 

lower triangularized stiffness matrix (figure 3) 

lower triangularized stiffness matrix of original structure (figure 3) 

load vector 

modified loading vector 
vector of reactions 
modified vector of reactions 

number of increments of stiffness changes (equation (39) ) 

upper triangularized stiffness matrix (figyre 3) 

displacement vector 

modified displacement vector 

multiplication factors (equation (39) ) 

increment 

shear stress 
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Subscripts 

av 

c 

eff 

f 

max 

r 

s 

u 


average 

modified degrees of freedom 
effective 

free degrees of freedom 
maximum 

modified constrained degrees of freedom 
constrained degrees of freedom 
unmodified degrees of freedom 


Superscripts 

(q) indicator of the number of stiffness increments (figure 3) 

T transpose of a matrix 


THE DIRECT MODIFICATION PROCEDURE 
The Structural Equations 


The usual governing equation for the static analysis of structures by the displacement 
method may be written as 


[K] {U} = {P} (1) 

where [K] denotes the structural stiffness matrix. {U} the unknown displacement vector and 
{P} the applied load vector. 

Equation (1), after elimination of the dependent equations (MPC) and including the forces 
of constraints, Qs, may be written in partitioned form as 
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where the subscripts 

f denotes the unconstrained (free) set 
and s denotes the constrained set. 

Expanding the simultaneous set of equations (2) and assuming that the enforced displacement 


vector {Ug} = 0 gives us 


1%J {Uf} = (Pf) 

(3) 

{Qs> = -{Ps> + [KftjT{Uf} 

(4) 






Equation (3) is solved for the unknown displacement vector {Uf } and then the forces of 
constraints {Qg}, are computed using equation (4). 

Preliminary Matrix Operations 


We define that 

c is the part of the f — set which is changed (modified) 
and u is the part of the f — set which remains unmodified. 


Equation (3) may be written in partitioned form as 
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Now if [Fff] is defined as the inverse of [K^f] , then the following relation is valid 
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Expanding equation (6) gives 


Kucl’^l^ucl + IKccl IFccI = tC.) 


*CC^ 


(7) 

( 8 ) 


By pre- and post-multiplying equation (7) by and respectively, the following 

is obtained 

[KuuJ'MKuJ = - [F,„] [F,,J-1 


uc-* cc 

By post-multiplying equation (8) by [F.,]‘^ and substituting equation (9) gives 


IFccl'^ = tK,,] 




uu' 


“uc^ 


(9) 


( 10 ) 


The effect of element modifications will be transmitted into [KffJ through a symmetric 
Boolean transformation matrix. The part of [Kff] affected by element modifications is limited 
to the symmetric submatrix [K^^] . Loading changes on the modified part of the structure are 
considered by the procedure. However, loading changes for the unmodified part of the structure 
should be considered through a usual restart job. 
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We now define the Boolean transformation matrices [3^,^] and [Byf]. Where 
[Bj,f] is the Boolean transformation matrix relating 

[Kff] to [K,c] 

[Buf] is the Boolean transformation matrix relating 
[Kff] to 1K„„1 

Each of the rows in [B^f] and [B^^'J contain all zeros, except for one unit value in the column 
relating to the position of the corresponding row in which is being changed or remains 


unmodified, respectively. 

Due to basic algebraic considerations the following orthonormal equations are valid 

[BcfJ [Bcf]T = (11) 

[BufJ IBufj’’ = [luu] (12) 

[Bcf^'^IBcfl + tBufl'^tBufl = [IffJ (13) 

The inverse of [Kc(] can now be related, via the Boolean transformation matrices, to its 

symmetric submatrices [F^^], [F^^,] and as follows 

[Fuul = IBuf) [Kffl-llBuflT ,14) 

IFucI = IBufl [KffJ-‘lBef]’' <I5) 

IFccl = [B(,fl lKffl-MBjf|T (16) 

The displacement vector {U(} is related, via the Boolean transformation matrices, to its 
subvectors {Upland {U^,} as follows 

{Uf} = [Buf]T{Uu> + [B^f]T{Uj;} • (17) 
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The Solution Procedure 


The equilibrium equations of the modified structures can be written in an equivalent form 
to equation (5) as 



where the element modifications in are expressed as 

IKccl = [K,,] + [A K,,] (19) 


[A Kj,gl is the stiffness matrix of the changes and may be singular. 


The final loading vector, {P^ }, on the modified part is expressed as 

{^c> = {Pc> + {^Pc> (20) 

where {A } defines the loading changes on the modified part of the structure. 


The final displacement vectors of the modified degrees of freedom, {U^,} ,and the 
unmodified degrees of freedom, {U^}, are expressed as 


(Uc> 

— 
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(21) 

{Uu> 
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(22) 


Where {A U,,} and {A U^} are the additional displacement vectors due to element modifications 
and loading changes. 

Using equations (19) to (22), and substituting for , {U^} , {0^} and {P^} in 
equation (18) and then subtracting equation (5) from it, yields 

I" 

[ AP, - AK^^ U, 

Using the first of equation (23) we get 

{AUu> = - [K^^] ^ {aU(,} (24) 
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Using the second of equation (23) and eliminating {AU^ } yields 

([AK^J + [K^cl - [Kuci'^[Kuur^[KuJ){Uc}= {APc>- IAK^cHUc) (25) 
By substituting equation (9) in equation (24) and equation (10) in equation (25) we obtain 

{iU„) = (F„J [K,(,1-'{4U,) (26) 

<[ IFcc'"') (('Uc> ' (6P^,}- (AK^J {Uj} (27) 

Using equation (17), the equivalent relation for the modified displacement vector 
{Uf + AUf} is as follows 

{Uf + AUf} = {Uu+AUu>+ [Bcfl'^ {Uc+AU^} (28) 

where {AUf} is the change of the displacement vector in the f-set. 


Subtracting equation (17) from equation (28) yields 

{AUf} = [Byf]T{AUu) + [Bcfl'^CAU^} (29) 

Substituting {AUy} from equation (26) in equation (29) yields 

(AUf) = IBufl’'[F„,)(F,jl-l (AU^) t [Bcf|T(AU,) (30) 

Using equation (15) and substituting for [F^g] in equation (30) we obtain 

{AUf} = [BufjT[Buf)[Kffl-l[BcflT[F^J{AUc} +[B^f]T{AUc} (31) 

Using equation (13) and substituting for ( [Puf^^f®uf^ ^ equation (31), yields 

{AUf} = ([Iff] - [Bcf]'^[Bcfl ) [Kff]-l[Bcf][Fccl-l {Uc}+[B^flT {AU^} 

(32) 

From equation (16) we establish the identity 

aBcf][Kff]-l[Bcf]T)[FcJ-l = [Iccl (33) 

Using the identity of equation (33) we reduce equation (32) to obtain 

{AUf} = [Kff]-l[Bcf]'^ [Fccl'^{AUc} (34) 
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Using equation (27) and substituting for {AU^,} in equation (34) we obtain the final 
equation for the change in the displacement vector of the f-set due to the modifications as 
follows 

{iUf) = ((K„| t (F„|-‘)-» / 

({AP,) - [AK^,1 {U,,}) Oi) 

And the modified displacement vector, {Uf } is expressed as | 

(Uf) = {Uf} + {AUf} (^5) 

Using equation (4) we define the equivalent relations for the modified vector of 
constraints, {Qg}, as 

{Qs> = - {Pg} + [iCfg]'^ {Uf} ( 37 ) 

where the modified matrix [K(g] is expressed as 

IKf,!’’ = lKf,|T + [B„)T[iK„lT|B^f] ,38) 

where 


lAKj.j.J is the compacted change matrix of [Kf^J and [B„J is the Boolean 
transformation matrix relating [Kf^] to [ AK^rJ , i.e., r denotes the part of the s - set 
which is changed. 

Considering the final equation for the modified displacement vector, {Uf } , (equation 36) 
we note that for the analysis of the modified structure, only [AK„] and {Ap’} have to be 

computed based on the modifications on the structure. However, the major computations are due 
to the following steps 

Performing a forward sweep and backward substitution, using the decomposed 
triangular matrices of [Kff] on the Boolean matrix [B^.f]T. 

Decomposition and performing a forward sweep and backward substitution using 

[Fee] on a unit matrix [leel . We note that [F ] is a positive definite symmetric 
matrix. 


The remaining computations are involved with smaller computing time. 
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The efficiency of the procedure is further improved if, after computing 
and (first two terms of equation (35)), changes are made only to the c degrees of freedom. 

This is implemented in the procedure described here by increasing or decreasing the term [AK^^l 
(last term of equation (35) ) by a constant factor as shown below: 

If 

ttj is an initial multiplication factor for [AK^^] and [AK^i-l 

a 2 is an incremental multiplication factor 

q is the number of incremental stiffness changes 

Then 

a=ttj+ q-a 2 

Therefore if [AK^j.] is multiplied within the computational loop by the variable a , then only 
the term ( [Fccl"^ recomputed for each increments of stiffness 

change. The factors a i , a 2 and q are defined in the NASTRAN analysis of the modified 
structure by standard PARAM bulk data cards. 

The generation of the Boolean transformation matrices [B^.^] and (Bj^l is tedious and 
error -prone. Therefore, a computer program, REAN, has been developed which generates those 
matrices automatically in a form suitable for NASTRAN. 

An overall schematic diagram of the procedure is shown in figure 1. 

DESCRIPTION OF PROGRAM REAN 

The REAN (RE ANalysis) computer program is a preprocessor for NASTRAN. It performs 
the following main tasks 

a) Updates MPT (Material Property Table) and EST (Element Summary Table) tables 

b) Generates the Boolean transformation matrices [B(,fj and [Bj.gl 

c) Generates the appropriate DMAP ALTER package for Rigid Format 1. 
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The execution of program REAM, illustrated in figure 2, is subdivided into the following steps. 

* The input to the program consists of 

/ 

(a) a standard NASTRAN Data Deck defining the grid points, element connections;' 
element properties, material properties for the structural modifications only and 
the values and number of incremental stiffness changes 

(b) the NASTRAN tables CSTM, YS, GPL, USET, OQGl, MPT and EST. 

* The program then checks 

(a) if CBAR, CQDPLT, CQUAD2, CTRIA2 or CTRPLT elements are to be modified; 
if so, the EST tables are updated. This is done in order to obtain the correct 
stresses for the elements listed above. 

(b) it the material properties defined for the modified part of the structure are 
different from those defined for the original structure; if so, the MPT tables 
are updated. 

Using the grid points defined for the modified part ot the structure in conjunction 
with the NASTRAN tables GPL, USET and OQGl the program first checks if any 
of the modified elements are connected to degrees of freedom which are constrained. 
If this is the case, the Boolean transformation matrices [B^j-] and [BJ. 5 I are 
generated. If none ot the modified elements are connected to constrained degrees of 
freedom, then only [B^f] is generated; [ B^f ] and [ Bj.^ ] are inserted via DM I 
card images. 

* Using the information obtained via the NASTRAN files UTland UT2 , the program 
then generates the appropriate DMAP Alter package and creates a NASTRAN input 
file. 

After the execution of program REAN, NASTRAN is loaded and the analysis for the modified 
structure is performed. 


DMAP ALTER PACKAGES 

The modification procedure is effected in NASTRAN utilizing standard DMAP statements 
The following two types of DMAP Alter Packages are required; 

1) The hrst DMAP Alter package, inserted in the original deck, utilizes OlTPUTl and 
OUTPUT 2 modules only. 
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OUTPUTl is used to save all the matrices and data blocks required to perform the 
analysis of the modified structure and the data recovery of nodal displacements, 
reactions, element forces and stresses. The OUTPUT2 module is used to save data 
blocks required for program REAN as explained in the previous section. 

2) The second DMAP ALTER package, generated by program REAN and which effects 
the solution procedure, is automatically inserted into the data deck of the modified 
structure and performs the actual reanalysis procedure. 

A flow diagram of the DMAP functional steps and their equivalent results is shown 
in figure 3. The DMAP statements for the first and second packages are listed and 
described in the Appendix. It should be noted that, according to the problem type, 
program REAN automatically chooses the appropriate form of inclusion of the 
optional data blocks CSTM, GM, KFSO, KSS, PS and YS in the SDRl and SDR2 
modules. 


NUMERICAL EXAMPLES 

The direct modification procedure has been applied to several real aircraft structures. 

Each of the examples described below is governed by different structural characteristics influ- 
encing the efficiency of the procedure. In order to achieve the minimum decomposition time for 
the regular analysis, resequencing of the grid points was performed. The examples were run on a 
CDC 6500 computer with 50000 jQ words of memory available for the execution of NASTRAN. 

Wing Structure (model 1) 

The wing structure was idealized using shear panels to represent the shear stiffness of the 
skins and webs, and rod elements representing their extensional stiffnesses. The mathematical 
model, figure 4, consisted of 6 BAR, 501 ROD, 288 SHEAR elements and 220 GRID points 
representing 1320 degrees of freedom of which 648 were unconstrained. The structure was sub- 
jected to two loading conditions and its stiffness matrix had 66 RMS column lengths. 

Three structural modifications were investigated; (a) the stiffness properties of 40 
elements were changed affecting 39 degrees of freedom; i.e. 6% of the structure; (b) the stiffness 
properties of 68 elements were changed affecting 85 degrees of freedom; i.e. 13% of the structure; 
and (c) the stiffness properties of 130 element were changed affecting 136 degrees of freedom; 
i.e. 21% of the structure. The changes were referenced as Regions “A”, “B” and “C”, respect- 
ively, on a plan view of the wing structure shown schematically in table 1 . This table also 
presents comparisons of the CPU execution time for a regular analysis versus analysis by the 
modification procedure. 


Wing Structure (model 2) 


This wing structure was a highly idealized mathematical model, as shown in figure 5. The 
structural elements used were the same as those used in model 1. The mathematical model con- 
sisted of 26 BAR, 1200 ROD, 660 CSHEAR elements and 560 GRID points representing 3360 
degrees of freedom, of which 1680 were unconstrained. The structure was subjected to two 
loading conditions and its stiffness matrix had 92 RMS column lengths. 

Three structural modifications were investigated: (a) The stiffness properties of 60 elements 
were changed affecting 66 degrees of freedom, i.e. 4% of the structure; (b) the stiffness proper- 
ties of 105 elements were changed affecting 132 degrees of freedom, i.e. 8% of the structure; 
and (c) the stiffness properties of 222 elements were changed affecting 23 1 degrees of freedom, 

I.e. 14% of the structure. The changes were referenced as regions “A”, “B” and “C”, respectively, 
on a plan view of the wing structure shown schematically in table 2. This table also presents com- ’ 
parison of the CPU execution time for a regular analysis versus analysis by the modification pro- 
cedure. 


Wing Structure (model 3) 

This model is similar in size to the one described as model 2, as shown in figure 5. How- 
ever instead of shear panels, isoparametric membrane elements, QDMEMl, were utilized. The 
mathematical model consisted of 42 CBAR elements, 662 QDMEMl, 860 ROD, and 
500 GRID points representing 3000 degrees of freedom, of which 1580 were unconstrained. The 
structure was subjected to three loading conditions and its stiffness matrix had 1 14 RMS column 
lengths. 


Three structural modifications were investigated: (a) The stiffness property of one element 
was changed affecting 6 degrees of freedom, i.e. 0.4% of the structure; (b) The stiffness proper- 
ties of 22 elements were changed affecting 48 degrees of freedom, i.e. 3% of the structure; and 
(c) the stiffness properties of 130 elements were changed affecting 225 degrees of freedom,' i.e. 
14% of the structure. The changes are referenced as Regions “A”, “B” and “C”, respectively, on 
a plan view of the wing structure shown schematically in table 3. This table also presents com- 
parisons of the CPU execution time for a regular analysis versus analysis by the modification pro- 
cedure. 


Fuselage Structure 


The fuselage was idealized using shear panels to represent the skins, rod elements to represent 
the stringers and beam elements to represent the bulkheads. The mathematical model, (figure 6) 
consisted of 160 BAR, 590 ROD, 350 SHEAR elements and 404 GRID points representing 2424 
degrees of freedom, of which 1293 were unconstrained. The structure was subjected to one loading 
condition and its stiffness matrix had 70 RMS column lengths. 


32 


Three structural modifications were investigated: (a) 3 elements were added to the 
structure affecting 28 degrees of freedom, i.e. 2.8% of the structure; (b) the properties of 12 
elements were changed affecting 72 degrees of freedom, i.e. 6% of the structure; and (c) the 
properties of 90 elements were changed affecting 180 degrees of freedom i.e. 14% of the struc- 
ture. The changes are referenced as Regions “A”, “B” and “C”, respectively, on an axonometric 
view of the fuselage structure shown schematically in table 4. This table also presents comparisons 
of the CPU execution time for a regular analysis versus analysis by the direct modification pro- 
cedure. 


Panel In Shear 

Figure 7 shows the effect of reduced thickness at the center of a rectangular panel on the 
modulus of rigidity and on the maximum shear stress. The panel was idealized by 25 rectangular 
QDMEM elements and the thickness of the center element was reduced to zero in 10 increments. 

This example demonstrates the procedure’s capability of “one-shot” solution of several 
cases differing in stiffness (i.e. with varying [AKj,j,] in equation (35) ). In addition, the re- 
quired CPU execution time utilizing the present procedure was 1/4 of the time required 
for 10 regular runs. 


DISCUSSION 

The efficiency studies for the solution times of the numerical examples solved, indicated 
the following intrinsic features of the direct modification procedure utilized. 

(a) For large problems whose stiffness matrices have large RMS values for their semi- 

bandwidth, the forward sweep and back substitution using the triangularized stiffness 
matrix of [K^f] on (First two terms of equation (35) ) is the most time- 

consuming. Table 3 shows that for 3% modification, for example, the solution time 
for the modification procedure is equal to the solution time for the regular analysis. 
However, for a structure with lower RMS values, the break-even point in the solution 
time is reached at approximately 5% modifications, as can be seen from table 1 . 

(b) For modifications affecting a relatively large number of degrees of freedom, the 

decomposition and subsequent forward sweep and back substitution using [Fj,^] on 
the unit matrix (third term of equation (35) ) may take a considerable portion 

of the solution time needed for the modification procedure. For example in table 1 , 
for modification “C”, the computation of [F^^,]"^ took 1/3 of the total solution 
time. 

(c) The types of elements used in the finite element model may have a considerable effect 
on the efficiency of the procedure. Table 3 shows that when changing the stiffness 
property of one element only 1/20 of the total CPU execution time of the regular 
analysis was required to obtain the new displacement vector. However the same 
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change for the same mesh with different types of elements would require 1/12 of the 
execution time of the regular analysis, as can be deduced from table 2. 

(d) When using the option of incremental stiffness changes, the efficiency of the procedure 
is increased substantially. The comparison between the regular analysis and the mod- 
ification procedure (table 5) shows that for all the cases solved, the break-even CPU 
time for solution of the displacement vectors came to more than 30% modifications 
of degrees of freedom. For Wing model 1 the break-even point was approximately 

60% modifications of degrees of freedom when including 4 additional stiffness incre- 
ments. 


CONCLUSIONS 

A direct modification procedure has been presented and its implementation in NASTRAN 
has been described. The efficiency of the procedure has been investigated by solving several real 
and large structures. The study of the relative CPU execution times for the various problems 
indicated that the efficiency of the procedure is dependent on the following main structural char- 
acteristics; (a) Size of the original structure and the number of active columns in the stiffness 
matrix; (b) the number of stiffness modifications; (c) the types of finite elements used in the 
mathematical model; and (d) the number of required additional stiffness modifications. 


The break - even point of the procedure varies between 8% to 60% modifications of the 
degrees of freedoms. However, the higher percentage can only be achieved when using the option 
of incremental stiffness modifications. It is therefore concluded that the direct modification 
procedure is highly efficient when investigating regions with continuous varying stiffness magnitudes, 
such as occurs during the post - buckling phenomena. Furthermore, the procedure was found to be' 
very useful when studying the effect of removing or adding elements between existing grid points 
to a structure or changing the stiffness properties of elements at different locations. 
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APPENDIX 


LISTING AND DESCRIPTION OF THE DMAP STATEMENTS 


DMAP Alter Statements For The Original Structure 


1. 

ALTER 

119 

2. 

OUTPUT 1 

LLL,ULV,KFS„//C,N,- 1/C,N,8/C,N,USERTPM $ 

3. 

OUTPUT 1 

USET,EQEXIN,SIL,BGPDT,//C,N,0/C,N,8 $ 

4. 

OUTPUT 1 

GM,CSTM,YS,PS,KSS//C,N,0/C,N,8 $ 

5. 

OUTPUT2 

MPT,EST,„//C,N,- 1/C,N,1 1/C,N,USERTP2 $ 

6. 

OUTPUT2 

CSTM,YS,GPL,USET,0QG1//C,N,12/C,N,USERTP3 

7. 

ENDALTER 



Description of DMAP Alter Statements 

2. The lower triangularized stiffness matrix of [Kj-j-J , the displacement vector and [Kfg] , 
which are used in the modification procedure, are written on file INP8. 

3, 4. NASTRAN tables to be used for the recovery of the nodal displacements, nodal reactions, 

element forces and stresses are written on file INP8. 

5. The Element Summary Table (EST) and the Material Property Table (MPT) are written on 
file UTl such that, if necessary, they can be updated by Program RE AN. 

6. Writes NASTRAN Tables on file UT2, with the help of which Program REAN computes 
the Boolean transformation matrices and generates the DMAP ALTER package for the 
subsequent analysis of the modified structure. 


DMAP Alter Package For The Analysis Of The Modified Structure 


1. ALTER 

2. INPUTTl 

3. INPUTTl 

4. INPUTTl 

5. INPUTT2 

6. FILE 

7. PARAM 


78 

/LLLO,ULVO,KFSO„/C,N,- 1 /C,N,8/C,N,USERTP1 $ 
/USETl .EQEXINl ,SIL1,BGPDT1 ,/C,N,0/C,N,8 $ 

/GM 1 ,CSTM 1 , YS 1 ,PS 1 ,KSS 1 //C,N,0/C,N,8 $ 

/MPTl ,EST1 ,„/C,N,- 1 /C,N, 1 1 /C,N,USERTP2 $ 

KFF 1 =SAVE/KFF2=SAVE/Q 1 =S AVE/XB=SAVE/V=SAVE/ 
BBT=SAVE/BS=SAVE/KFSA= SAVE /KFSB=SAVE $ 

//C,N,NOP/V,N,NC=-l $ 
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8. 

PARAM 

//C,N,NOP/V,N,NA=-l $ 

9. 

PARAM 

//C,N,SUB/V, Y,NCASES/V,Y,NCASES= I /C,N, 1 $ 

10. 

LABEL 

W3 $ 

11. 

COND 

W7,NC $ 

12. 

COND 

W1,NA $ 

13. 

LABEL 

W7 $ 

14. 

PARAM 

//C,N,ADD/V,N,NC/C,N,1 $ 

15. 

FBS 

LLLO„BBT/V/C,N, 1 /C,N, 1 /C,N, 1 /C,N, 1 $ 

16. 

MPYAD 

BB,V,/Q/C,N,0/C,N, 1 /C,N,0/C,N, 1 $ 

17. 

MPYAD 

BB,BBT,/DD/C,N,0/C,N, 1 /C,N,0/C,N, 1 $ 

18. 

MPYAD 

BB,ULVO,/XB/C,N,0/C,N, 1 /C,N,0/C,N, 1 $ 

19. 

DECOMP 

Q/QL,QU/C,N,1/C,N,1 $ 

20. 

FBS 

QL„DD/QI/C,N, I /C,N, 1 /C,N, 1 /C,N, 1 $ 

21. 

ADD 

KFS,/KFSA/C,Y,FACA = (0.0, 0.0) $ 

22. 

ADD 

KFS,/KFSA/C,Y,FACB = ( 1 .0,0.0) $ 

23. 

ADD 

KFF,/KFF1/C,Y,FACA = (O.O.O.O) $ 

24. 

ADD 

KFF,/KFF 1 /C, Y,FACB = ( 1 .0,0.0) $ 

5 **** TOP OF DMAP LOOP FOR INCREMENTS OF STIFFNESS CHANGES **** 

25. 

LABEL 

W1 $ 

26. 

ADD 

KFSA,KFSB/KFSC $ 

27. 

EQUIV 

KFSC,KFSA/NA $ 

28. 

ADD 

KFF1,KFF2/KFF3 $ 

29. 

EQUIV 

KFF3,KFF1/NA $ 

30. 

PARAM 

//C,N,ADD/V,N,NC/V,N,NC/C,N,1 $ 

31. 

PRTPARM 

//C,N,0/C,N,NC $ 

32. 

MPYAD 

KFF3,XB,/XBB/C,N,0/C,N, 1 /C,N,0/C,N, 1 $ 

33. 

ADD 

QI,KFF3/R $ 

34. 

DECOMP 

R/RL,RU/C,N,I/C,N,I $ 

35. 

FBS 

RL„XBB/U1/C,N,1/C,N,1/C,N,1/C,N,1 $ i 

36. 

SMPYAD 

V,QI,U1,„/U2/C,N,3/C,N,-I/C,N,1 $ I 

37. 

ADD 

U2,ULVO/ULVN $ 

38. 

JUMP 

LBL9 $ 

39. 

ALTER 

110,119 

40. 

SMPYAD 

BBT,KFSC,BS,„KFSO/KFSl/C,N,3/C,N,l/C,N,l $ 
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1 

1 

I 

1 

1 


41. 


SDRl 


USET 1 „ULV 1 „YS 1 „GM 1 ,PS1 ,KFS 1 ,KSS 1 ,/UGV,PGG,QG/ 


42. 

43. 

44. 

45. 
$ 

46. 


V,N,NSKIP/C,N, STATICS $ 

SDR2 CASECC,CSTMM,MPT1„EQEXIN1,SIL1„,BGPDT1,PGG,QG, 

UGV.ESTl ,/OPGl ,0QG1,0ES1 ,0EF1 ,PUGV1 /C,N, STATICS $ 

ALTER 121 

PARAM //C,N,SUB/V,Y,NCASES/V,Y,NCASES/C,N,1 $ 

COND W3,15 $ 

**** BOTTOM OF DMAP LOOP FOR INCREMENTS OF STIFFNESS CHANGES 
ENDALTER 




Description Of DMAP ALTER Statements 

1. The procedure starts after, the stiffness matrix for the modified structure, [AK^^js] has 

been assembled. 

2. Retrieves the lower triangularized stiffness Matrix of [K^f] , LLLO, Displacement vector, 
ULVO, and KFSO for the unmodified structure. 

3. 4. Retrieves NASTRAN tables to be used for the recovery of the nodal displacements, reactions, 

element forces and stresses. It should be noted that not all the data block names saved via 
OUTPUTl in the original analysis appear here. Program REAN checks if the files saved are 
empty, retrieving only those data blocks which contain information required for the analysis 
of the modified structure. 

5. Retrieves the updated Material Property Table (MPTl) and Element Stiffness Table (ESTl). 
The update of the tables is performed by REAN. 

7—9. Defines the default values for a ^ , a 2 and q. 

15. Performs a forward sweep and backward substitution using the lower triangular matrix of 

[Kff] , on this yields 

16. Multiplies [B^.f] • • [B^f]'^ which yields the matrix [F^.^] • 

18. Multiplies [B(,fl by the displacement vector ULVO of the unmodified structure. This yields 
the displacement vector, {U^.} , for the part of the structure which is to be modified. 

19. Decompose the matrix into an upper and lower triangular matrix. 

20. Performs a forward sweep and backward substitution using the lower triangularized matrix 
of [Fgj.] on a unit matrix, obtained by multiplying [B^f] • [B^f]^, this yields [F^.^] 
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21-24. 


25. 

32. 

33. 

34. 

35. 

36. 

37. 

40. 

41, 42. 


Initializes [AK^^] and [AK^j.] of the modified structure. 


Start of the DMAP loop for the specified increments of stiffness changes. 
Multiplies [AK^^j • {U^.} 

Adds ( [AK^J + ) 

Decomposes the matrix ( [AK,,^.] + ) into upper and lower triangular 

matrices. 


Performs a forward sweep and backward substitution using the lower triangular matrix 
of ( [AK^J + [F^^]-l ) on ( [AK^,J {U^} ) 

This yields ( [AK^J + [F^J-1)-1 • [aK^^] • {U^} 

Multiplies the results obtained in 15, 20 and 35. This yields the displacement vector, 

{AUi"}, which represent the influence of the structural modification on the original 
structure. 

Adds ULVO and {AUf}and yields the final displacement vector ULVN of the modified 
structure. 

Computes the nodal reaction (See equation (37) ). The Boolean transformation matrix 
is computed by REAN. 

DMAP Statement 40 is not included by Program REAN if none of the modified elements 
are incident at degrees of freedom which are constrained. 

Recovers via SDRl and SDR2 the nodal displacements, reactions element stresses and 
forces for the modified structure. 


45. 


Bottom of DMAP loop. 
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Table 1 —Comparison of CPU Execution Times for Wing Model 






































Table 3 —Comparison of CPU Execution Timjes for Wing Model 













Table 4 -Comparison of CPU Execution Times for Fuselage 
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FIGURE 1: OVERALL SCHEMATIC DIAGRAM FOR 1 CYCLE REANALYSIS 
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FIGURE 3: FLOW DIAGRAM FOR BLOCK 2 
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FIGURE 4: IDEALIZATION OF AIRCRAFT WING (MODEL 1) 
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FIGURE 6: IDEALIZATION OF AIRCRAFT WING (MODELS 2 and 3) 





